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Abstract 

We calculate the current and differential conductance for the junction between a 
superconducting (SC) STM tip and a Luttinger liquid (LL). For an infinite single-channel 
LL, the SC coherence peaks are preserved in the tunneling conductance for interactions 
weaker than a critical value, while for strong interactions {g < 0.38), they disappear and 
are replaced by cusp-like features. For a finite-size wire in contact with non-interacting 
leads, we find however that the peaks are restored even for extremely strong interactions. 
In the presence of a source-drain voltage the peaks/cusps split, and the split is equal 
to the voltage. At zero temperature, even very strong interactions do not smear the 
two peaks into a broader one; this implies that the recent experiments of Y.-F. Chen 
et. al. (Phys. Rev. Lett. 102, 036804 (2009)) do not rule out the existence of strong 
interactions in carbon nanotubes. 

Scanning tunneling spectroscopy (STM) is becoming an important tool for accessing the 
electronic properties of low- dimensional systems. Thus, in the past, the local density of states 
(DOS) for various materials such as the cuprates graphene |2], and carbon nanotubesIH] 
has been studied, both theoretically and experimentally. In general the STM tip is assumed 
to be a non-interacting metal with a constant DOS, and the system to be analyzed is in 
equilibrium (the voltage is constant throughout the sample). This allows one to extract the 
unknown DOS directly from the STM tunneling conductance. 

Recently a few studies have also concentrated on tips that are not normal metals, and 
that may have an intrinsic variation in the DOS, such as superconducting (SC) STM tips 
in [5]. These experiments have the potential of measuring not only the DOS, but also other 
quantities such as the Fermi distribution. The first such experiment P| looked at a non- 
interacting system, and found that in equilibrium the differential conductance from a SC 
STM tip shows the characteristic SC gap and coherence peaks, while in the presence of a 
voltage bias each peak splits into two peaks, with the distance between them being given 
by the applied voltage. For a non-interacting system this splitting can be traced back to an 
out-of-equilibrium double-step Fermi distribution. 

Similar experiments have been performed recently on carbon nanotubes [5j, but the in- 
terpretation of these measurements is not so straightforward, due to the presence of strong 
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electronic interactions which make the elementary excitations no longer fermionic but rather 
fractionally-charged. One might naively expect that the Fermi distribution is ill-defined in 
these systems, and that the strong electronic interactions contribute to the smearing of the 
two SC coherence peaks into a single one. On the other hand the experiments of |5| show the 
presence of two coherence peaks, which, when combined with this naive expectation, appear 
to imply that the interactions in nanotubes are weak ! 

To fully assess the implications of this experiment, and to test whether the expectation 
that the peaks merge is borne out of rigorous calculations, one needs to study theoretically the 
injection of electrons into a strongly-interacting out-of-equilibrium one-dimensional system. 
This is rather a challenging problem, that has not begun to be addressed theoretically until 
quite recently P [3 [HI [9] . 

We study the injection from a SC tip into an out-of-equilibrium LL, and we show that, 
even in the presence of very strong interactions, the tunneling conductance is very similar to 
that of a non-interacting wire. Our calculations therefore demonstrate that the apparently 
non-interacting features observed in |5j do not rule out the existence of strong interactions in 
carbon nanotubes. 

We begin by analyzing an infinite single-channel LL wire in equilibrium, and we find that 
the tunneling conductance displays SC coherence peaks for interactions weaker than a critical 
value, corresponding to a fractional charge parameter, g, larger than Qc = 0.38. Beyond this 
value the coherence peaks disappear, and are replaced by cusp-like features. Upon taking into 
account the finite size of the wire [101 HH [12] and the presence of non-interacting contacts (as 
in the experiment of [5]) we find that the peaks reappear for any value of g. Furthermore, 
the tunneling conductance shows also oscillations with a period inverse proportional to the 
length of the wire and has a background power-law dependence whose exponent depends on 
the strength of the interactions. 

If a voltage is applied between the ends of the wire, the peak/cusp features split. The 
magnitude of the split is exactly equal to the applied voltage. At zero temperature we never 
see a smearing of the two peaks into a single one, regardless of the strength of the interaction. 
Hence, the expectation that interactions cause the two peaks to merge into one is not correct^. 
This in turn implies that the presence of two peaks found in [S] can not be used to rule out 
the presence of strong interactions in nanotubes ! 

Our calculations show that the experimental setup of [5J can be used to asses the strength 
of the interactions if one focuses instead on the power-law dependence of the density of states. 
If one measured the density of states for a wide-enough range of voltages (up to tens of meV) 
the data should allow to extract the value of the interaction parameter g. 

A quantum wire connected to metallic leads is described by the Hamiltonian 

n = no + nv. (i) 

Here Tio describes the interacting wire and the leads in the framework of the inhomogeneous 
single-channel LL model [10], |1T[ [121 U\ (we will discuss at the end how our results are affected 

^This is consistent with the theoretical calculations in Ref. [8J showing a double-step Fermi distribution in 
an out-of-equilibrium LL. 
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by the presence of the extra channels of conduction in nanotubes). Exphcitly: 
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(2) 
(3) 



where Tiv describes the chemical potential applied to the wire. The interaction parameter g{x) 
is space-dependent and its value is g in the bulk of the wire, and 1 in the leads. For convenience, 
the end-points of the wire are denoted by Xi = —L/2 and X2 = L/2. The chemical potential 
is chosen such that fi{x) = —Vsd/'^ for x < —L/2 (left lead), fj,{x) = Vsd/'^ for x > L/2 (right 
lead), and /i(a;) = for -L/2 < x < L/2 

Tunneling is allowed between the wire and a superconducting tip at at x = 0. A schematic 
view of the system is shown in Fig. [31 The voltage of the tip is fixed at V. 
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Figure 1: A quantum wire adiabatically coupled to metallic leads. The leads are held at different 
chemical potentials fii = -eVso/'^ and ^2 = ^ysn/'^- The electrons can tunnel from and into a 
superconducting STM tip at x = 0. The voltage of the tip is fixed at V . 



The Hamiltonian for the SC tip is assumed to be of the BCS type with a linear dispersion 
(efc oc k) and with a SC BCS coupling. 



Htip = X] X] ekc\^Cka + (Acfc|C_fc| + h.c). 

k a=T,i 



(4) 



Our results do not depend on the details of the dispersion and of the dimensionality of the 
tip, as long as the tip DOS has the standard BCS form z/t(e) = 9{e — A)e/V — A^. 
The tunneling between the tip and the wire is described by the local Hamiltonian 



Ht = tciix = 0)[^f (x = 0) + ^^(x = 0)] + h.c 



(5) 



where \I/q''^(x) are the chiral fermionic interacting operators that can be related to the free 
bosonic modes of the LL via bosonization [13] . 
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We use the non-equilibrium Keldysh formalism and we perform a perturbative expansion 
up to second order in the tunneling coefficient t, using the formalism developed in Ref. [7J to 
calculate the tunneling current between a normal tip and an out-of-equilibrium LL. We find 
that the tunneling current between the SC tip and the wire at zero temperature is 

POO 

I oc / dtQTit)lm{exp[n{t) + iI{t)]J^sit)} , (6) 
Jo 

where for a finite-size wire jTj 

VsDt 



Qrit) = 2 cos sm{Vt) , (7) 



and 
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with r = tvp/gL, p = (1 — g)/{l + g), and aw is a dimensionless short-length cutoff. The 
quantity Qrit) incorporates the effects of the applied tip voltage and source-drain voltage, 
while TZit) and I{t) are related to the Green's functions of the system and incorporate infor- 
mation about the finite size via r. For an infinite wire one should substitute these quantities 
by: 

nt) = -^{g + g-')^n^^^ (14) 



167r afy 



^(^) = -^(9 + g'' - 2) arctan ( — ) . 



(15) 



Here aw is a short-time (high-energy) cutoff for the infinite wire. 

For the superconducting tip, by taking a Fourier transform of the usual SC Green's function 
we obtain 

J^sit) = AK,[{tt + aT)A], (16) 
where qt is a short-time (high-energy) cutoff for the tip, and Ki is the first Bessel K function. 
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For a finite wire we can use the formulae above to compute tlie differential conductance 
dl/dV numerically. For an infinite wire this calculation can be performed analytically, by 
rewriting (along the lines of [8]) the derivative of the tunneling current as: 



dl ^ . r , dufie 

r]=±l 



V |tp r de^^ile + V + vVsd/2\'^ -\e~V- vVsnm (17) 
Jo 



where 'y = {g+l/g — 2)/2 and z/((e) is the density of states of the superconducting tip. This 
integral now gives: 

^ oc tJ2 l^ni^ vVsD/2r'/\F^i3/2, -1/2-, 3/2 - r, a/2 + V/2 + vVsd/^) + 

\A + V + vVsD/2r'^'2F,{3/2, -1/2; 3/2 - 7; A/2 - V/2 - vVsd/2)] (18) 

where 2-^1 is the corresponding hypergeometric function. 

We expand this function (in equilibrium - Vsd = 0) close to the singularity at V = A and 
at V > A. When ^ A: 

§^^{V- A)-V2 . (19) 

For non-interacting systems 7 = 0, and the tunneling conductance diverges as expected: 
^ (X (y — A)~^/^. In the presence of interactions this exponent is modified. For very strong 
interactions 7 increases, the exponent changes sign, and the divergence is replaced by a power- 
law-type cusp. This regime is achieved when 7 > 1/2, which corresponds to g < 0.38. 
When \^ ^ A we obtain 

§^V^. (20) 

For non-interacting systems the tunneling conductance is constant at large V. In the presence 
of interactions the conductance acquires a power-law background with a positive exponent 7 
as expected for tunneling into a LL. 

In Figs. 2 a) and b) we plot the tunneling conductance for Vsd = 0, for g = 1 and g = 0.5. 
We take A = 1 in arbitrary units. In Figs. 2 c) and d) we plot the corresponding tunneling 
conductanc^ when Vsd = A. As described above, for g = 0.5 > gc the SC coherence peaks 
are preserved in the spectrum. When the wire is out of equilibrium, Eq. (IT5]) implies that the 
tunneling conductance is given by the sum of two tunneling conductances, which are shifted 
respectively by ±Vsd/2. 

In Figs. 3 a) and c) we repeat the same calculation for g = 0.3 < gc- We see that indeed the 
peaks disappear and are replaced by cusp-like features, corresponding to a power-law decay 
of the tunneling conductance ffTI?]) . 

Our results change significantly if we take into account the finite size of the wire and the 
presence of metallic contacts. The integrals can now only be done numerically. In Fig. 3 b) 
and d) we plot the tunneling conductance for g = 0.3 for a finite-size wire with characteristic 
finite-size energy scale = hvp/gL = A/3. We see that the peaks that disappear at infinite 
lengths are restored for a finite-size system. This is consistent with previous observations for 



^The source-drain voltage is applied such that the average of the right and left potentials is zero, and the 
bulk of the wire remains neutral. If the voltage of the bulk of the wire becomes non-zero (because of, say, a 
gate voltage), the curves will shift along the voltage axis such that the symmetry point of the curve will no 
longer be at zero tip voltage, but at a value of the tip voltage equal to the gate voltage. 
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Figure 2: The tunneling conductance from an STM tip into an infinite Luttinger liquid (in arbitrary 
units) for LL interaction parameters g = 1 (in a) and c)) and g = 0.5 (in b) and d)). The source-drain 
voltage VsD is zero for a) and b), and is V^d = A = 1 for c) and d). 



nanotubes flUi [TTl [T2] : when V — A the physics is dominated by long wavelengths and 

is hence non-interacting. 

Most importantly, for both infinite and finite systems, the two features at A ± ^5/5/2 are 
not smeared into a single one, regardless of the strength of the interactions. Thus, as far as the 
peaks are concerned, even a very strongly interacting wire behaves as non-interacting. This is 
consistent with previous observations of a double-step out-of-equilibrium Fermi distribution 
[8], and explains why the features observed in Ref. [5] appear for both interacting and non- 
interacting systems. The only signatures of interactions are the disappearance of the peaks at 
very small g for the infinite wire, and the power-law dependance of the tunneling conductance. 

We should note that for a nanotube the parameter 7 is reduced from {g + 1/g — 2)/2 to 
[g + 1/ g — 2)/8 because of the extra three channels of conduction [T^. Hence, the critical 
value of g below which the peaks disappear is ^^cn ~ 0.17. Nevertheless, the finite-size effects 
are always relevant for non-chiral LLs, and hence we expect that the spectra presented in 
Figs. 3 b) and d), characterized by oscillations and SC coherence peaks, should describe the 
tunneling conductance between a SC tip and a nanotube in all situations. 

For a physical set of parameters, A = ^ = 1.5meV and hvp/gL ^ 0.5meV (corresponding 
to a tube length of ~ 4/im at g = 0.25), our results for the tunneling conductance reproduce 
very well those measured in [5] (compare for example Figs. 3 b) and 3 d) here with Figs. 2 c) 
and 3 in [5]). 

To summarize, we have calculated the dependence on voltage of the STM tunneling con- 
ductance between a superconducting STM tip and an out-of-equilibrium Luttinger liquid. We 
have found that for an infinite single-channel LL in equilibrium, one can observe SC coher- 
ence peaks for interactions weaker than a critical value [g < 0.38). Beyond this value, the 
coherence peaks disappear and are replaced by cusp-like features. We have also found that 
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Figure 3: The tunneling conductance from a SC STM tip into an infinite (a) and c)) and finite-sized 
(b) and d)) g = 0.3 Luttinger liquid. The source-drain potential Vsd is zero for a) and b) and is 
VsD = A = 1 for c) and d). The energy corresponding to the finite size of the wire, hvp/gL, is taken 
to be A/3. 



the peaks are restored if one takes into account the finite size of the LL and the presence of 
non-interacting contacts. For both infinite and finite-sized wires, the tunneling conductance 
shows a background power-law dependence, with an exponent that can be used to determine 
the strength of the interactions. 

In the presence of an applied voltage between the ends of the wire, the tunneling-conductance 
features (peaks or cusps) split in two, and the magnitude of the split is equal to the voltage. 
Thus, there is no smearing of two peaks into a single one, regardless of the strength of in- 
teractions. This implies that the measurements of Ref. [5], that reveal features which appear 
naively to be non-interacting, do not in fact rule out the presence of strong interactions in 
carbon nanotubes. 
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